topological invariants in two dimensions from quantum Monte Carlo 
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We employ quantum Monte Carlo techniques to calculate the Z2 topological invariant in a two- 
dimensional model of interacting electrons that exhibits a quantum spin Hall topological insulator 
phase. In particular, we consider the parity invariant for inversion-symmetric systems, which can 
be obtained from the bulk's imaginary-time Green's function after an appropriate continuation 
to zero frequency. This topological invariant is used here in order to study the trivial-band to 
topological-insulator transitions in an interacting system with spin-orbit coupling and an explicit 
bond dimerization. We discuss the accessibility and behavior of this topological invariant within 
quantum Monte Carlo simulations. 
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I. INTRODUCTION 

Topological insulators have been intensively explored 
in recent yearsi*^, especially since their prediction^ and 
experimental realization in HgTe quantum wells^. Pro- 
posed also by Kane and Mele in a theoretical model of 
spin-orbit interactions in graphene^ in search for an in- 
trinsic quantum spin Hall (QSH) effect, topological band 
insulators in two dimensions are characterized in the 
presence of time-reversal invariance by a Z2 topologi- 
cal index of the insulating electronic statQ^^. In the case 
of non-interacting systems^, the Z2 topological invariant 
can be extracted from the insulating band structure in 
analogy to the Thouless, Kohmoto, Nightingale and den 
Nijs (TKNN) classification of Block wave functions rel- 
evant for the integer quantum Hall effeclii. Still in the 
context of non-interacting systems, it was found, that 
for inversion symmetric systems, e.g., in the sublattice- 
symmetric case on the graphene lattice, the Z2 topo- 
logical invariant can be easily extracted directly from 
the Hamiltonian matrix of the system at the so-called 
time-reversal invariant momenta (TRIM) in the Brillouin 
TjOm^—. At these specific momenta, Kramers degener- 
ate partners share the same band-structure eigenvalues, 
and from the parity of the occupied band eigenstates the 
corresponding Z2 parity invariant (PI) is easily obtained. 
This approach will be reviewed below within a more gen- 
eral setting. 

Recently, topological insulators augmented with 
(strong) electron-electron interactions have attracted 
growing attention (see e.g. Ref. [ll| for a recent review 
of work on two-dimensional systems). Hence, the ques- 
tion arises, how the concept of a topological characteri- 
zation of an insulating electronic state can be extended 
beyond the non-interacting band-structure regime. An 
important issue is, how such topological information can 
be efficiently calculated for interacting systems, in par- 
ticular using unbiased numerical methods, such as quan- 
tum Monte Carlo (QMC) simulations. Several means to 



calculate topological invariants for interacting electronic 
insulators have been put forwardiSrJi. In a non-trivial 
generalization from the non-interacting case, these topo- 
logical quantities are constructed based on the system's 
dressed single-particle Green's function, which remains 
a well-defined quantity also for interacting systems. Of 
particular interest from a numerical perspective are the 
schemes presented in Refs. [l5| and[l^, which allow one to 
obtain the topological index based solely on the single- 
particle Green's function G'(w, k) at zero frequency cj = 
and momentum k. As will be shown below, this quan- 
tity can be easily obtained from QMC calculations. The 
calculations can be further simplified for systems with 
explicit inversion symmetry, where G{ui = 0, k) needs 
to be obtained at the TRIM onlyi^, similar to the non- 
interacting case^. A description of this approach to ex- 
tract the PI for interacting systems, which in addition 
also exhibit spin S*^ conservation, will be presented be- 
low. 

Such methods to obtain topological invariants for in- 
teracting systems have been applied recently, e.g., to cor- 
related electron systems in one dimension using the nu- 
merically exact time-dependent density matrix renormal- 
ization group (DMRG) approachi^. For two-dimensional 
interacting systems, approximate means to estimate the 
Green's function have been employed; for example the 
variational cluster approximation (VGA) has been ap- 
plied to the Kane-Mele model with local interactions (the 
Kane-Mele-Hubbard modeli^i"— )^ to study the transi- 
tion from the quantum spin Hall topological insulating 
phase to the antiferromagnetic Mott insulator regime 
at strong interactions. Dynamical mean-field theory 
(DMFT) has been employed to study the interaction- 
driven transition between topological states in a Kondo 
insulator— and cluster DMFT to study the three dimen- 
sional pyrochlore iridates^i. 

Here, we set out to employ unbiased and numerically 
exact methods to access Green's function-based topolog- 
ical invariants in two dimensional fermion systems. In 
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FIG. 1. (Color online) (a) Honeycomb lattice with a unit 
cell indicated by dashed lines. The arrows indicate the lat- 
tice vectors ai,2. Bold (red) lines denote bonds with hopping 
amplitude t' , while the hopping amplitude along the other 
bonds on the honeycomb lattice equals t. Filled (open) cir- 
cles indicate lattice sites belonging to the A (B) sublattice. 
(b) Brillouin zone with the time-reversal invariant momenta 
{r, Ml, M2, M3}, and the reciprocal lattice vectors bi,2 indi- 
cated. 



particular, we use a projective QMC scheme to study 
the PI for a Kane-Mele-Hubbard model with anisotropic 
hopping, which exhibits a topological insulator regime, a 
trivial, non-magnetic insulating phase, as well as an an- 
tiferromagnetically ordered Mott insulating regime. We 
analyze the PI in these phases, the transitions between 
them and assess the Pi's characterization of these differ- 
ent regimes. The goal of this paper is not to provide a 
detailed analysis of the complete phase diagram of this 
model, but to instead illustrate the actual application of 
the Green's function approach to study topological in- 
variants in interacting two-dimensional fermion systems. 

The rest of this paper is organized as follows: In 
the next section, we introduce the dimerized Kane- 
Mele-Hubbard model that we explore further below, fol- 
lowed by a review of how to extract the PI from the 
zero-frequency Green's function. We examine the non- 
interacting limit of our model, where the Green's func- 
tion and the PI may be easily calculated. After that, we 
discuss how to obtain the PI for finite interactions from 
QMC simulations, before we then apply this approach to 
the dimerized Kane-Mele-Hubbard model. 



the spin-orbit ncxt-nearest-neighbor term 

and the Hubbard local interaction term 

Hu = U ^{a\^ai^a\^a,^ + b\^b,^b\^bi^) . 

i 

Here, (6-^.) denote creation operators for spin-cr 
fermions on a sublattice- A site (sublattice-B site), with i 
denoting the two-site unit cell at position on the honey- 
comb lattice. The spin-orbit coupling strength is denoted 
by A, while Uij = ±1 depending on whether the consid- 
ered hopping process involves a left, or a right turn. We 
allow for a different nearest-neighbor hopping strength t' 
along one of the three nearest-neighbor bond directions, 
as compared to the other directions, cf. Fig. 1. Here, 
the unit cell is chosen such that it contains a t' bond and 
is centered on this bond. The two lattice vectors of the 
honeycomb lattice = ao(3/2, ±-\/3/2) are also shown 
in Fig. 1(a). In the following, we set the distance between 
nearest neighboring lattice sites ag = 1. For t' = t, the 
usual Kane-Mele-Hubbard model is recovered, which for 
finite spin-orbit coupling A and in the small- C/ regime 
features a QSH topological insulating region, adiabati- 
cally connected to the U = Q QSH state. Increasing the 
onsite repulsion U eventually drives the system into an 
ordered phase with long-ranged transverse antiferromag- 
netic correlations^!. Furthermore, at [/ = 0, the explicit 
bond dimerization allows to drive the system from the 
topological insulator QSH state to a (trivial) band in- 
sulating phase for t' > 2t. At t' /t = 2, the system is 
gapless with the bulk gap closing at one of the M-points 
in the Brillouin zone, cf. Fig. 1(b). This will be exam- 
ined in more detail below as well as the properties of the 
model for t' > t and finite interactions, [/ > 0. To study 
the effects of interactions in terms of the topological in- 
variants, we employ quantum Monte Carlo simulations 
to calculate the imaginary-time Green's functions of this 
model Hamiltonian and then transform to the Green's 
function at zero frequency, from which we extract the Z2 
PI for this inversion symmetric system. 



II. DIMERIZED KANE-MELE-HUBBARD 
MODEL 

In the following, we consider the half-filled Kane-Mele- 
Hubbard modeli^"— 1^ with an additional, explicit bond 
dimerization, described by the Hamiltonian 

H = Ho + Hso + Hu, (1) 
with the nearest-neighbor hopping terms 

i a (ij) a 



III. PARITY INVARIANT FROM GREEN'S 
FUNCTION 

In an inversion symmetric system, the PI may be 
calculated from the system's Green's function follow- 
ing Ref. [isl which generalizes the procedure from the 
non-interacting casG^. Here, due to the explicit S*^ con- 
servation of the Hamiltonian, the Green's function is 
block-diagonal in spin-space, and the procedure can be 
restricted to a single spin sector. The zero-frequency 
Green's function for each spin sector, 6*^(0, k), where 
a = +1 (—1) for spin-up (spin-down), thus is a 2 x 2 
matrix in the A/B-sublattice basis. Denoting by bi^2 
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the reciprocal lattice vectors (with b; • a.j = 27r5jj), we 
consider the four TRIM 

Kni,n2 ="lbl/2 + ?l2b2/2, 71, = 0, 1 , (2) 

corresponding to the F-point and the three M-points in- 
dicated in Fig. 1(b), at which the operation of inver- 
sion commutes with the zero-frequency Green's function 
Go- (0,k). Here, the operation of inversion that inter- 
changes the two sublattices and squares to the identity 
can be represented in the sublattice basis by the first 
Pauli- matrix, P ~ a^. Simultaneously diagonalizing the 
two matrices P and G(j(0, /«„j^_„2), we identify for each of 
the four TRIM the eigenvalue of P for the common eigen- 
vector with a positive eigenvalue of Go-(0, H;ni,n2 

). These 

eigenvectors are referred to as right-zeros or R-zeros in 
Ref. |T^. Denoting the corresponding P eigenvalue of the 
R-zero by t]k.„^ „, , we obtain the PI, A e {0, 1}, as 

(-1)^ - n . (3) 

ni ,712 

from any of the two spin sectors, which together form a 
Kramer's pair at each TRIM. It is thus sufficient for the 
calculation of the PI, to only consider e.g. the spin-up 
sector due to the explicit conservation of the Hamil- 
tonian. The procedure is however easily generalized to 
inversion symmetric systems without explicit S*^ conser- 
vationi^. 



IV. NON-INTERACTING CASE 

To illustrate the above procedure, let us first consider 
the non-interacting limit, i.e., the dimerized Kane-Mele 
model. For U = 0, the Hamiltonian H can be diagonal- 
ized directly via a transformation to momentum space, 

^ = E(4,.0'^^(k) (foH- (4) 

In each spin sector a = +1 ( — 1), the Hamiltonian matrix 
at wave vector k equals 

where gw = t' + t{e''^^'^ + e''*^ '') relates to the near- 
est neighbor hopping terms and 7k = 2A(— sin(\/3fcy) + 
2cos(3fca./2) sin(A/3fc.y/2)) to the spin-orbit term. The 
system described by H conserves Sz, such that the 
Green's function G(aj, k) is block-diagonal in spin-space, 
and each spin component in the non-interacting case 
equals 

G,{uj,k) = [uj - h,{-k)r\ (6) 

At zero frequency this is essentially the inverse of the 
Hamiltonian matrix: 

GAO,k) = -h-\k). (7) 
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FIG. 2. (Color online) Band structure of the dimerized Kane- 
Mele model along the indicated path through the Brillouin 
zone for A/t = 0.2 and different values of t', as indicated. 



Based on the approach outlined in the previous section, 
we then obtain for finite values of A a change in the PI 
from A = 1 for i' < 2t to A = for t' > 2t. This indi- 
cates the change from a topological insulator to a trivial 
band insulating state driven by the explicit bond dimer- 
ization. At t' = 2t, the system becomes semi-metallic 
due to the single particle gap closing at the M3 point, 
i.e., at k = Ki^i. This can be seen from the band struc- 
ture shown for X/t = 0.2 in Fig. 2. In the following, we 
will examine this transition also at finite values of U. 
Before performing such an analysis, we first explain, how 
we extract the PI in the interacting regime from QMC 
simulations. 



V. PARITY INVARIANT FROM QMC 

Once the zero-frequency Green's function G(0,k) has 
been obtained for the interacting model, the PI can 
be calculated as outlined in Sec. IIIII In analogy with 
Eq. ([7]) for the non-interacting case, one can associate 
to the interacting model a fictitious Hamiltonian ma- 
trix ft,topoi(k) = — G~^(0,k), which has been dubbed the 
topological HamiltoniariF^. It contains the topological in- 
formation of the interacting model, where for the free case 
/Jtopoi(k) equals the Hamiltonian matrix of H. Hence, we 
merely need to consider, how the zero-frequency Green's 
function is obtained from the QMC calculations. In par- 
ticular, we employed a projective QMC scheme, by which 
we obtain the momentum and spin resolved single par- 
ticle Green's function in imaginary time within the sys- 
tem's ground state on finite lattices. To obtain the PI, 
we then calculate from the imaginary-time data of the 
Green's function those at Matsubara frequencies, and 
continue in particular to zero frequency. For this pur- 
pose, let us first consider the system at a finite temper- 
ature T = 1//3. The imaginary time Green's function 
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Gcr(T, k;/3) at a given momentum k and spin projection 
cr is a two-by-two matrix with entries 



[G.(r,k;/3)],, --(ck,.,,(T)ct (0));3, 



(8) 



where j,l = 1,2 is a sublattice index, with Ck,(T,i = o-k.a 
and Ck,cr.2 = &k,(T- For frequencies a;„ = 2{n + 1)tt/ /S the 
Matsubara-Green's function is then given as 



G,(iw„,k;/3) = / Ga(r,k;/3)e'""^dr. 



(9) 



Particle hole symmetry of the model at half- filling, i.e., 

under the transformation cj^.^^- — > (ik,crj = (^l)"'cLk o- j 
in each spin sector together with inversion symmetry 
leads to the following conditions on G^r (r, k; /3) : For equal 
sublattices, [Go-(r, k; = [Ga{l3 — t, —k; /3)]jj . while, 
for j^/, [GAr,k-p)y=-[G,{p-T,-k;l3)],i. 

We thus obtain for the diagonal elements of the Green's 
fimction at one of the TRIM n = Kni,n2 the equation 



[Ga (iw„ , k; 13)] jj = 2 i / [Gair^n; /3)] ^ j sin(a;„T) dr , 
Jo 

(10) 

and, for j ^ I, 

[Gai'iujn, K,; /3)]j7 = 2 / [Gct(t, k; cos(w„t) dr . 
Jo 

(11) 

Now, the limit (3 ^ oo can be taken properly: From 
the projective QMC, we obtain the ground state Green's 
function Ga{T,K,) = lim^_j.oo Gct(t, k; /?), and then per- 
form the above integrals with a sufficiently large cutoff 
/3 — )■ 0, set e.g. by the imaginary time evolution length 
of the Green's function employed in the QMC simu- 
lations. Here, we used 6 = 20/i. This cutoff proved 
to be sufficient for the Greens function to decay to zero 
within error bars, especially for large values of U /t, but 
for the extreme cases close to the topological-to-trivial 
band insulator transition, where the gap becomes very 
small. Note, that one cannot simply take the limit 
ibjn — before accounting for the (anti) symmetry con- 
ditions on the imaginary time Green's functions. This 
would lead to wrong results, as exemplified below. Af- 
ter (anti)symmetrization, the limit iw„ — ^ can be per- 
formed with the T = Green's functions, so that in 
particular. 



and, for j ^ Z, 



[G^{lo = Q,k)]j, =0, 



61/2 



(12) 



[G,(c^ = 0,k)],, =2 / [G,{T,K)]jidT. (13) 
Jo 

Hence, within the QMC simulations, one merely needs 
to measure the off-diagonal part of the Green's function 
explicitly. To illustrate the above point, consider for a 



moment the non-interacting limit, for which the exact 
T = imaginary-time Green's function 



G,{r,K) = ~\e-\^-\-(\ ^ 



(14) 



If calculated naively, via Go-(t, F) dr, one would 
(wrongly) obtain a finite value of [Gct(w = 0,F)]jj in- 
stead of the actual value (i.e. zero), which also follows in 
this case directly from Eq. ([6]). 



VI. QMC RESULTS 

After having examined the calculation of the PI for 
the interacting system in the previous section, we now 
present results from QMC simulations of the dimer- 
ized Kane-Mele-Hubbard model. We employ a projector 
axillary-field determinantal QMC scheme^^ by which we 
obtain the momentum and spin resolved single particle 
Green's function in imaginary time within the system's 
ground state for finite lattices with N = 2L^ lattice sites 
employing periodic boundary conditions. Here, L de- 
notes the linear system size, which for multiples of six al- 
lows all TRIM as well as the corners of the Brillouin zone 
(the so-called Dirac points) to be presented. In particu- 
lar, we use a projection length Q — 50/t, imaginary-time 
step At = 0.05/t and linear systems sizes L = 6, 12 and 
18. An imaginary time evolution length = 20/t has 
been used to obtain the Green's function, as discussed in 
Sec. V. Details on the employed QMC method in appli- 
cation to the Kane-Mele-Hubbard model can be found in 
Ref. m. 

To test the feasibility of extracting the PI within QMC, 
we first consider the t'/t-driven transition between the 
topological insulator regime and the trivial band insula- 
tor for large t' at finite values of U . In the following, 
we consider \/t = 0.2, in order to focus on the QSH 
to dimerized insulator transition without being compro- 
mised by the influence of the QSH-insulator transition 
at A = 0, and without loss of generality<2ii2i22ii2£ As an 
example. Fig. [3] shows the imaginary-time dependence of 
the off-diagonal component of the Green's function at the 
Ms-point at Kh, which in the following we denote by 



Go(r) [Gt(T,Kn) 



(15) 



Indeed, a change in the PI in our model can be traced 
back to a sign-change in Go(t) [more precisely, in the 
corresponding integral of Eq. (|13p ]. As can be seen from 
Fig. 131 for U /t ~ 2 and \/t = 0.2, this change occurs 
between t' /t ~ 1.94 and t' /t = 1.96, and correspond- 
ingly, A jumps from A 1 to A = between these 
values. This indicates, that for these parameters, the 
topologieal-to-trivial band insulator transition occurs for 
a slightly smaller values of t' /t = 1.95(1) than at [/ 0, 
where the transition takes place at precisely t' /t — 2. 
This can be understood to be the consequence of the 
super-exchange induced by the local Coulomb repulsion 
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FIG. 3. (Color online) Off-diagonal component of the Green's 
function at the Ma-point for L = &, U /t = 2 and X/t = 0.2 at 
various values of t'/t. Error bars are of the order of the line 
width and have been omitted for clarity. 
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FIG. 4. (Color online) Evolution of the single particle gap 
Aap for different system sizes as a function of t'/t for U/t = 
2 and X/t = 0.2 near the quantum phase transition from 
the topological insulator to the dimerized phase. The inset 
focuses in on the transition region. 

which favors the singlet formation on the i'-bonds. At 
the transition point, the single-particle excitation gap 
Asp closes, as can be seen from Fig. 21 which shows 
Asp at the M-point Kh, obtained from the decay in 
imaginary time of the diagonal Green's function elements 
[G-|-(r, Kii)]jj cx cxp (— rAgp). This reflects the same gap 
closing at the transition point as observed for [/ = at 
t'/t = 2. 

While in the small- 1/ region, the QSH state is stable 
with respect to interactions and adiabatically connected 
to the U = limit, for sufficiently large values of U the 
system enters a transverse antiferromagnetically ordered 
Mott-insulating phase, where the time-reversal symme- 
try of the Hamiltonian H is spontaneously broken in the 
thermodynamic limiti^Ti^ii^iSl. This transition is how- 



ever not related to a closing of the single particle gap, as 
has been demonstrated by unbiased QMC simulations. 
The single particle gap only exhibits a local minimum at 
the transition point, but does not close2ii2^. This result 
from numerically exact simulations is in contrast to pre- 
vious VGA calculations, which concluded that the single 
particle gap closes at the transition to the antiferromag- 
netic phased. In fact, the Green's function exhibits no 
qualitative change across the transition. This can be seen 
also from the QMG data in Fig. [SJ where Go (t) is shown 
for different values oi U a.t X/t = 0.2 and for t' — t. 

From previous QMG simulations2i, we know that long- 
ranged antiferromagnetic order sets in for these param- 
eters near U/t ~ 5 and flux induced edge states are 
absent^. However, Go{t) exhibits no significant changes 
in this interaction region. In particular, and in contrast 
to the t'-scan considered above, Go(t) does not exhibit a 
change in its sign. That this is not a finite size effect, can 
be seen in the inset of Fig. [SJ where we compare QMC 
data at C//t = 8 for two different system sizes, L = 6 and 
L = 12, which are seen to indeed be finite-size converged. 
We verified that also up to U/t = 40, no sign change oc- 
curs in Go (t) . This implies that the PI A stays constant 
when tuning across the antiferromagnetic transition. We 
verified explicitly, that even at A = the PI takes on a 
non-trivial value in the antiferromagnetic Mott insulating 
region. 

How does this relate to the quantum phase transition 
that takes place when the system enters the antiferromag- 
netic region, which is thus not adiabatically connected to 
the U = state? Only in the thermodynamic limit anti- 
ferromagnetic order persists, which spontaneously breaks 
time-reversal and the inversion (sublattice) symmetry of 
the Hamiltonian. Yet this is not monitored by the single 
particle Green's function, on which the calculation of the 
PI is based. Spontaneous symmetry breaking in the or- 
dered region implies a degenerate ground state subspace 
in the thermodynamic limit. In each specific ground state 
from this manifold, the sublattices A and B arc not equiv- 
alent anymore, and this condition for a well defined PI 
is broken. Remarkably, even in the antiferromagnetic re- 
gion, the degeneracy of the ground state manifold implies 
the existence of low-energy gapless excitations, namely 
the Goldstone modes. However, these soft spin excita- 
tions are of particle-hole type, and thus not attainable in 
the single-particle sector. 

Interestingly, once the antiferromagnetic order is de- 
stroyed in the system by a sufficiently strong dimeriza- 
tion t', the PI does change to a trivial value. To ex- 
plore this behavior, let us start from the large-t' region, 
t' > 2t. Starting at C/ = from the trivial band in- 
sulator region for t' > 2t and switching on local inter- 
actions U > 0, the system remains insulating, and also 
does not develop long-range order. This can be most 
directly seen in the large-JJ limit. Here, the effective 
model for the low-energy physics is a Heisenberg model 
with an exchange dimerization along the <'-bonds: The 
exchange interaction J' = At''^ /U in second order pertur- 



6 



o 



0.50 



0.40 



0.30 - 



0.20 



0.10 



0.00 




o 




FIG. 5. (Color online) Off-diagonal component of the Green's 
function at the Ms-point for L — 6, \/t = 0.2, t' = t and 
different values of U from U = to U = 8t (top to bottom). 
Across the transition from the QSH insulator to the antiferro- 
magnetic insulator the Green's function remains qualitatively 
unchanged. Error bars are of the order of the line width and 
have been omitted for clarity. Inset: The Green's function 
shows very little finite size dependence due to the large gap 
induced by the large coupling parameters. 



FIG. 6. (Color online) Off-diagonal component of the Green's 
function at the Ma-point for L = 6, U/t = 8 and X/t = 0.2 
at different values of t' between t' = t and t' = 1.6t (top to 
bottom). Error bars are of the order of the line width and 
have been omitted for clarity. 



as across the transition. 



bation theory is more than a factor of 4 larger than the 
exchange interaction J = At^ /U along the other nearest- 
neighbor bonds, and also dominates over the weak (for 
\/t — 0.2) next-nearest-neighbor anisotropic exchange 
coupling J2 = 4A^/C/ related to the spin-orbit tcrmsi^. 
The strong J'-dimerization drives the spin system into a 
non-magnetically ordered, dimerized phase. Indeed, for 
J2 = 0, the critical ratio beyond which the antiferro- 
magnetic order vanishes in the Heisenberg model on the 
dimerized honeycomb lattice equals J'/J = 1.735(1)^, 
which relates here to a ratio of t' /t w 1.32 in the Hub- 
bard model in the large-L/ limit. For t' > 2t, the system 
thus resides inside a non-magnetic phase, adiabatically 
connected to the trivial band insulator at [/ = 0. Corre- 
spondingly, the PI of the system does not change upon 
increasing U at fixed t' > 2t. 

On the other hand, decreasing the ratio t' /t at suffi- 
ciently large U, a transition from the large-t' non-ordered 
phase to the antiferromagnetic phase occurs, and we ob- 
serve a corresponding transition in the PI: In Fig. [51 we 
consider in particular the case oi U/t = 8 and X/t = 0.2. 
Upon varying t' /t, we find a change in the sign of Go{t), 
and in more detail, the PI changes from A = 1 to A = 
beyond t' /t = 1.28(2). Remarkable is the fact, that this 
transition is again related to the emergence of low-energy 
gapless excitations, namely the Goldstone modes, which 
appear in the antiferromagnetic phase, but not in the 
dimerized phase, where instead a finite spin-gap sepa- 
rates the singlet ground state from the lowest triplet ex- 
cited state, that relates in the strong- J' limit to a triplet 
excitation on one of the strong J' bonds. The single par- 
ticle gap Asp however stays finite in both phases, as well 



VII. CONCLUSIONS 

We explored Green's function based methods to obtain 
topological invariants in a two-dimensional strongly in- 
teracting fermion system that exhibits trivial band insu- 
lating, Mott insulating and topological insulator regimes. 
Given an adiabatic connection for a phase of the inter- 
acting system to the non-interacting limit, we found that 
the calculated parity invariant indeed does not change, 
and thus allows extracting the direct transition between 
the topological insulator and the trivial insulator region 
for finite interactions. However, since the parity invari- 
ant relates to the single particle Green's function, and 
hence captures single particle properties only, it does not 
allow monitoring, e.g., the transition to the antiferro- 
magnetic regime from the topological insulator side. A 
change of the parity invariant would require correspond- 
ing changes in the single particle Green's function, which 
are not being observed in this case. Use of the parity 
invariant based on approximate methods to calculate the 
Green's function may however lead to deviating conclu- 
sions. For example, the change in the parity invariant 
within the variational cluster approximation to the Kane- 
Mcle-Hubbard model^^, is accompanied by a closing of 
the gap in the single particle Green's function^, which 
does not reflect the actual behavior of this mode P^'^'^ . 
Although its usage is thus restricted, we have shown that 
the parity invariant nevertheless constitutes a readily ac- 
cessible measure within quantum Monte Carlo simula- 
tions for a large variety of (quantum) phase transitions 
from topological to trivial insulators. 

Recently we became aware of a QMC investigation^ 
which examined the parity invariant in a related model, 
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focusing on the interaction region below the magnetic 
ordering transition. 
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